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This is a brief review of recent theoretical efforts to understand persistence in nonequilibrium 
systems. Some of the recent experimental results are also briefly mentioned. I also discuss recent 
generalizations of persistence in various directions and conclude with a summary of open questions. 



The problem of persistence in spatially extended 
nonequilibrium systems has recently generated a lot of in- 
terest both theoretically and experimentally [j|4l0| . 
Persistence is simply the probability that the fluctuating 
nonequilibrium field does not change sign upto time t. 
These systems include various models undergoing phase 
separation process [p]4i| Jll] , [l^ , [Tr^ , ^4] , p^| , simple diffusion 
equation with random initial conditions |5||| , several re- 
action diffusion systems in both pure |l7| and disordered 
|L8| environments, fluctuating interfaces [34H36I, Lotka- 



Volterra models of population dynamics 1 19 and granular 
medium pc[ . 

The precise definition of persistence is as follows. Let 
4>{x, t) be a nonequilibrium field fluctuating in space and 
time according to some dynamics. For example, it could 
represent the coarsening spin field in the Ising model after 
being quenched to low temperature from an initial high 
temperature. It could also be simply a diffusing field 
starting from random initial configuration or the height 
of a fluctuating interface. Persistence is simply the prob- 
ability Po(t) that at a fixed point in space, the quantity 
sgn[(f>(x,t) — ((f>(x,t))] does not change upto time t. In 
all the examples mentioned above this probability decays 
as a power law Po(t) ~ t~ e at late times, where the per- 
sistence exponent 9 is usually nontrivial. 

In this article, we review some recent theoretical efforts 
in calculating this nontrivial exponent in various models 
and also mention some recent experiments that measured 
this exponent. The plan of the paper is as follows. We 
first discuss the persistence in very simple single vari- 
able systems. This makes the ground for later study of 
persistence in more complex many body systems. Next 
we consider many body systems such as Ising model and 
discuss where the complexity is coming from. We follow 
it up with the calculation of this exponent for a sim- 
pler many body system namely diffusion equation and 
see that even in this simple case, the exponent 6 is non- 
trivial. Next we show that all these examples can be 
viewed within the general framework of the "zero cross- 
ing" problem of a Gaussian stationary process (GSP). We 
review the new results obtained for this general Gaussian 
problem in various special cases. Finally we mention the 
emerging new directions towards different generalizations 
of persistence. 

We start with a very simple system namely the one 



dimensional Brownian walker. Let </>(t) represent the po- 
sition of a 1-d Brownian walker at time t. This is a single 
body system in the sense that the field <fi has no x depen- 
dence but only t dependence. The position of the walker 
evolves as, 



(i) 



where 77(f) is a white noise with zero mean and delta cor- 
related, (<j)(i)(j)(t')) = S(t - t'). Then persistence P (t) 
is simply the probability that <f>(t) does not change sign 
upto time t, i.e, the walker does not cross the origin upto 
time t. This problem can be very easily solved exactly 
by writing down the corresponding Fokker-Planck equa- 
tion with an absorbing boundary condition at the origin 
pjj . The persistence decays as Po(t) ~ t -1 / 2 and hence 
9 = 1/2. The important point to note here is that the 
exact calculation is possible here due to the Markovian 
nature of the process in Eq. (1). Note that (j> evolves 
according to a first order equation in time, i.e., to know 
4>{t), we just need the value of <j>(t — At) but not on 
the previous history. This is precisely the definition of a 
Markov process. 

In order to make contact with the general framework 
to be developed in this article, we now solve the same 
process by a different method. We note from Eq. (1) 
that 77(f) is a Gaussian noise and Eq. (1) is linear in <f>. 
Hence (f> is also a Gaussian process with zero mean and a 
two time correlator, (cf>(t)(j>(t r )) = min(t,t') obtained by 
integrating Eq. (1). We recall that a Gaussian process 
can be completely characterized by just the two-time cor- 
relator. Any higher order correlator can be simply calcu- 
lated by using Wick's theorem. Since min(t, t') depends 
on both time f and t' and not just on their difference 
|t — t'|, clearly <f> is a Gaussian non-stationary process. 
From the technical point of view, stationary processes 
are often preferable to non-stationary processes. Fortu- 
nately there turns out to be a simple transformation by 
which one can convert this non-stationary process into 
a stationary one. It turns out that this transformation 
is more general and will work even for more complicated 
examples to follow. Therefore we illustrate it in detail for 
the Brownian walker problem in the following paragraph. 

The transformation works as follows. Consider first the 
normalized process, X(t) = </>(f)/ sj (<f> 2 (f)). Then X(t) is 
also a Gaussian process with zero mean and its two-time 
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correlator is given by, (X(t)X(t')) = min(t,t')/ \/(tt'). 
Now we define a new "time" variable, T = log(f). Then, 
in this new time variable T, the two-time correlator be- 
comes, (X(T)X(T')) = exp(-|T - T'|/2) and hence is 
stationary in T. Thus, the persistence problem reduces 
to calculating the probability Po(T) of no zero crossing 
of X(T), a GSP characterized by its two-time correlator, 
(X(T)X(T'))=eM-\T-T'\/2). 

One could, of course, ask the same question for an arbi- 
trary GSP with a given correlator (X(T)X(T')) = f(\T- 
T'\) [in case of Brownian motion, f(T) = exp(-T/2)]. 
This general zero crossing problem of a GSP has been 
studied by mathematicians for a long time 0] . Few re- 
sults are known exactly. For example, it is known that 
if f(T) < 1/T for large T, then P (T) - exp(-^T) 
for large T. Exact result is known only for Markov 
GSP which are characterized by purely exponential cor- 
relator, f(T) = exp(-AT). In that case, P (T) = 
■|sin _1 [exp(— AT)] [p2[ . Our example of Brownian mo- 
tion corresponds to the case when A = 1/2 and therefore 
the persistence Pq(T) ~ exp(— T/2) for large T. Revert- 
ing back to the original time using T = log(i) , we recover 
the result, Po(t) ~ t~ x / 2 . Thus the inverse of the decay 
rate in T becomes the power law exponent in t by virtue 
of this "log-time" transformation. Note that when the 
correlator f(T) is different from pure exponential, the 
process is non-Markovian and in that case no general an- 
swer is known. 

Having described the simplest one body Markov pro- 
cess, we now consider another one body process which 
however is non-Markovian. Let <p(t) (still independent of 
x) now represents the position of a particle undergoing 
random acceleration, 

J=V(t) (2) 

where r)(t) is a white noise as before. What is the prob- 
ability Po(£) that the particle does not cross zero upto 
time t? This problem was first proposed in the review 
article by Wang and Uhlenbcck |2^] way back in 1945 
and it got solved only very recently in 1992, first by Sinai 
0, followed by Burkhardt [gU by a different method. 
The answer is, Po(t) ~ t" 1 / 4 for large t and the persis- 
tence exponent is — 1/4. Thus even for this apparently 
simple looking problem, the calculation of is nontrivial. 
This nontriviality can be traced back to the fact that this 
process is non-Markovian. Note that Eq. (2) is a second 
order equation and to know <p(t + At) , we need to know 
its values at two previous points (f>(t) and 4>{t — At). Thus 
it depends on two previous steps as opposed to just the 
previous step as in Eq. (1). Hence it is a non-Markovian 
process. 

We notice that the Eq. (2) is still linear and hence </>(t) 
is still a Gaussian process with a non-stationary correla- 
tor. However, using the same T = log(t) transformation 
as defined in the previous paragraph, we can convert this 



to the zero crossing problem in time T of a GSP with 
correlator, f(T) = §exp(-T/2) - ±exp(-3T/2). Note 
that this is different from pure exponential and hence 
is non-Markovian. We also notice another important 
point: It is not correct to just consider the asymptotic 
form of f(T) ~ | exp(— T/2) and conclude that the ex- 
ponent is therefore 1/2. The fact that the exponent is 
exactly 1/4, reflects that the "no zero crossing" proba- 
bility Po(T) depends very crucially on the full functional 
form of f(T) and not just on its asymptotic form. This 
example thus illustrates the history dependence of the 
exponent 9 which makes its calculation nontrivial. 

Having discussed the single particle system, let us now 
turn to many body systems where the field <f)(x, t) now 
has x dependence also. The first example that were stud- 
ied is when cj)(x, t) represents the spin field of one dimen- 
sional Ising model undergoing zero temperature coarsen- 
ing dynamics, starting from a random high temperature 
configuration. Let us consider for simplicity a discrete 
lattice where <fr(i,t) = ±1 representing Ising spins. One 
starts from a random initial configuration of these spins. 
The zero temperature dynamics proceeds as follows: at 
every step, a spin is chosen at random and its value is up- 
dated to that of one of its neighbours chosen at random 
and then time is incremented by At and one keeps repeat- 
ing this process. Then persistence is simply the probabil- 
ity that a given spin (say at site i) does not flip upto time 
t. Even in one dimension, the calculation of Po(t) is quite 
nontrivial. Derrida et. al. [Q solved this problem exactly 
and found P (t) ~ for large t with = 3/8. They also 
generalized this to g-state Potts model in 1-d and found 
an exact formula, 0(g) = — g + -^[cos -1 -^ — q)/qy/2}} 2 
for all q. 

This calculation however can not be easily extended 
to d = 2 which is more relevant from experimental point 
of view. Early numerical results indicated that the ex- 
ponent ~ 0.22 H for d = 2 Ising model evolving with 
zero temperature spin flip dynamics. It was therefore 
important to have a theory in d = 2 which, if not exact, 
at least could give approximate results. We will discuss 
later about our efforts towards such an approximate the- 
ory of Ising model in higher dimensions. But before that 
let us try to understand the main difficulties that one 
encounters in general in many body systems. 

In a many body system, if one sits at a particular point 
x in space and monitors the local field <p(x,t) there as a 
function of t, how would this "effective" stochastic pro- 
cess (as a function of time only) look like? If one knows 
enough properties of this single site process as a function 
of time, then the next step is to ask what is the prob- 
ability that this stochastic process viewed from x as a 
function of t, does not change sign upto time t. So the 
general strategy involves two steps: first, one has to solve 
the underlying many body dynamics to find out what the 
"effective" single site process looks like and second, given 
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this single site process, what is its no zero crossing prob- 
ability. 

Before discussing the higher dimensional Ising model 
where both of these steps are quite hard, let us discuss a 
simple example (which however is quite abundant in na- 
ture) namely the diffusion equation. This is a many body 
system but at least the first step of the two-step strategy 
can be carried out exactly and quite simply. The second 
step can not be carried out exactly even for this sim- 
ple example, but one can obtain very good approximate 
results. 

Let (f>(x, t) (which depends on both x and t) denote 
field that is evolving via the simple diffusion equation, 
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(3) 



This equation is deterministic and the only randomness 
is in the initial condition <p(x,0) which can be chosen as 
a Gaussian random variable with zero mean. For exam- 
ple, <fi(x, t) could simply represent the density fluctuation, 
<fi(x, t) = p{x, t) — (p) of a diffusing gas. The persistence, 
as usual, is simply the probability that <p(x,t) at some 
x does not change sign upto time t. This classical diffu- 
sion equation is so simple that it came as a surprise to 
find that even in this case, the persistence Po(t) ~ t~ e 
numerically with nontrivial 9 w 0.1207, 0.1875, 0.2380 in 
d = 1, 2 and 3 respectively. 

In light of our previous discussion, it is however easy 
to see why one would expect nontrivial answer even in 
this simple case. Since the diffusion equation (3) is lin- 
ear, the field 4>{x, i) at a fixed point a; as a function of t is 
clearly a Gaussian process with zero mean and is simply 
given by the solution of Eq. (3), (f>(x,t) — J d d x'G(x — 
x«,t)(f>(x',0), where G{x,t) = {Ant)-^ 2 exp[-a; 2 /4t] is 
the Green's function in d. Note that by solving the 
Eq. (3), we have already reduced the many body dif- 
fusion problem to an "effective" single site Gaussian pro- 
cess in time t at fixed x. This therefore completes the 
first step of the two-step strategy mentioned earlier ex- 
actly. Now we turn to the second step, namely the 
"no zero crossing" probability of this single site Gaus- 
sian process. The two time correlator of this can be 
easily computed from above and turns out to be non- 
stationary as in the examples of equations (1) and (2). 
However by using the T = log(i) transformation as be- 
fore, the normalized field reduces to a GSP in time T 
with correlator, {X{Tx)X{T 2 )) = [sech(T/2)] rf / 2 , where 
T = \Ti — T%\. Thus once again, we are back to the 
general problem of the zero crossing of a GSP, this time 
with a correlator f(T) = [sech(T/2)f/ 2 which is very 
different from pure exponential form and hence is non 
Markovian. The persistence, Po(T) will still decay as 
P {T) ~ exp(-0T) - t-° for large T (since f(T) decays 
faster than \/T for large T) but clearly with a nontrivial 
exponent. 



Since persistence in all the examples that we have dis- 
cussed so far (except the Ising model) reduces to the zero 
crossing probability of a GSP with correlator f(T) [where 
f(T) of course varies from problem to problem], let us 
now discuss some general properties of such a process. 
It turns out that a lot of information can already be 
inferred by examining the short-time properties of the 
correlator f(T). In case of Brownian motion, we found 
f(T) = exp(-T/2) ~ 1 - T/2 + 0(T 2 ) for small T. For 
the random acceleration problem, f(T) — | exp(— T/2) — 
±exp(-3T/2) - 1 - 3T 2 /8 + 0(T 3 ) for small T and 

for the diffusion problem, f(T) = [sech(r/2)] d / 2 - 1 - 
d_ T 2 + 0(T 3) ag T o In general /(T) = 1 — aT a + ... 

for small T, where < a < 2 [^2). It turns out that 
processes for which a = 2 are "smooth" in the sense that 
the density of zero crossings p is finite, i.e., the number 
of zero crossings of the process in a given time T scales 
linearly with T. Indeed there exists an exact formula 
due to Rice p = ^-f"{Q)/n when a = 2. How- 
ever, for a < 2, /"(0) does not exist and this formula 
breaks down. It turns out that the density is infinite for 
a < 2 and once the process crosses zero, it immediately 
crosses many times and then makes a long excursion be- 
fore crossing the zero again. In other words, the zero's 
are not uniformly distributed over a given interval and 
in general the set of zeros has a fractal structure |^6| . 

Let us first consider "smooth" processes with a = 2 
such as random acceleration or the diffusion problem. It 
turns out that for such processes, one can make very good 
progress in calculating the persistence exponent 9. 

The first approach consists of using an "independent 
interval approximation" (IIA) Consider the "effec- 
tive" single site process 4*{T) as a function of the "log- 
time" T = \og{t). As a first step, one introduces the 
"clipped" variable a — sgn(cf)), which changes sign at the 
zeros of 4>(T). Given that 4>{T) is a Gaussian process, it 
is easy to compute the correlator, A(T) = (er(O)cr(T)) = 
f sin" 1 [/ (T)] , where / (T) is the correlator of 4>{T) . Since 
the "clipped" process <r(T) can take values ±1 only, one 
can express A(T) as, 



A(T) = j2(-i) n Pn(n 



(4) 



where P n {T) is the probability that the interval T con- 
tains n zeros of 4>(T). So far, there is no approximation. 
The strategy next is to use the following approximation, 

P n (T) = (Ty 1 ( <ffi f 2 dT 2 ... [ dT n 

Jo JTx JT„_! 

x Q(Ti)P(T 2 - Ti) . . . P(T„ - T n _!)Q(T - T n ), (5) 

where P(T) is the distribution of intervals between two 
successive zeros and Q(T) is the probability that an in- 
terval of size T to the right or left of a zero contains no 
further zeros. Clearly, P{T) = -Q'(T). (T) = l/p is 



3 



the mean interval size. We have made the IIA by writ- 
ing the joint distribution of n successive zero-crossing 
intervals as the product of the distribution of single in- 
tervals. The rest is straightforward ||. By taking the 
Laplace transform of the above equations, one finally ob- 
tains, P(s) = [2 - F(s)]/F(s) where, 



F(a) = l + ^-a[l-sA(a)], 
2p 



(6) 



where (n r ) c are the cumulants of the moments. Using 
p = 1 — e, we express the right hand side as a series 
in powers of e. Fortunately the computation of the mo- 
ments of n is relatively straightforward, though tedious 
for higher moments. We have already mentioned the re- 
sult of Rice for the first moment. The second moment 
(ri 2 ) was computed by Bendat ]28| ]. We have computed 
the third moment as well ||. For example, for 2-d diffu- 
sion equation, we get the series, 



where the Laplace transform A(s) of A(T) can be eas- 
ily computed knowing f(T). The expectation that the 
persistence, Pq(T) and hence the interval distribution, 
P(T) ~ exp(— 9T) for large T, suggests a simple pole in 
the P(s) at s = —9. The exponent 9 is therefore given 
by the first zero on the negative s axis of the function, 

F{s) = 1 + - - / dTexpi-sTjsm-^fiT)}}. (7) 

For the d iffusion equation, f(T) = [sech(T/2)] d / 2 
and p = \J d/Sn 2 . We then get the IIA estimates 
of 6 = 0.1203, 0.1862 and 0.2358 in d = 1, 2 and 
3 respectively which should be compared to the sim- 
ulation values, 0.1207 ± 0.0005, 0.1875 ± 0.0010 and 
0.2380 ± 0.0015. For the random acceleration problem, 
f(T) = |exp(-T/2) - ±exp(-3T/2) and p = V3/2tt 
and we get, 9u a = 0.2647 which can be compared to its 
exact value, 9 = 1/4. 

Note that the IIA approach, though it produces ex- 
cellent results when compared to numerical simulations, 
cannot however be systematically improved. For this 
purpose, we turn to the "series expansion" approach ||] 
which can be improved systematically order by order. 
The idea is to consider the generating function, 



p(j>,t) = j2 P n p n (t) 



(8) 
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where P n (t) is the probability of n zero crossings in time 
t of the "effective" single site process. For p = 0, P(0, t) 
is the usual persistence, decaying as t~ e ^ as usual. Note 
that we have used 8(0) instead of the usual notation 9, 
because it turns out || that for general p, P(p, t) ~ 
for large t, where 9(p) depends continuously on p for 
"smooth" Gaussian processes. This has been checked nu- 
merically as well as within IIA approach Q. Note that 
for p = 1, P(l,t) = 1 implying 9(1) = 0. For smooth 
Gaussian processes, one can then derive an exact series 
expansion of 9(p) near p = 1. Writing p n = exp(nlogp) 
and expanding the exponential, we then obtain an ex- 
pansion in terms of of moments of n, the number of zero 
crossings, 



log P(p,t) 



(9) 



9(p = l-e) = ±e+(±-±y + 0(e 3 ). 

ZTT 7T Z 47T 



(10) 



Keeping terms up to second order and putting e = 1 
(in the same spirit as e expansion in critical phenomena) 
gives, 6»(0) = (7r + 4)/47r 2 = 0.180899 . . ., just 3.5% below 
the simulation value, 9 S i m — 0.1875 ± 0.001. This thus 
gives us a systematic series expansion approach for calcu- 
lating the persistence exponent for any smooth Gaussian 
process. 

Note that both the above approaches (IIA and series 
expansion) are valid only for "smooth" Gaussian pro- 
cesses (a = 2) with finite density p of zero crossings. 
What about the nonsmooth processes where < a < 2, 
where such approaches fail? Even the Markov process, 
for which f(T) = exp(— AT) is a non-smooth process with 
a = 1. Fortunately however for the Markov case, one 
knows that the persistence exponent 9 = X exactly. One 
expects therefore that for Gaussian processes which may 
be nonsmooth but "close" to a Markov process, it may be 
possible to compute 9 by perturbing around the Markov 
result. 

In order to achieve this, we note that the persistence 
Po(T) in stationary time T, can be written formally 0] 
as the ratio of two path integrals, 



Po(T) 



2^ >o mr)eMS} _ Zi 
JV(/>(T)exp[-S\ Z 



(11) 



where Z\ denotes the total weight of all paths which never 
crossed zero, i.e., paths restricted to either positive or 
negative (which accounts for the factor 2) side of <fi = 
and Zq denotes the weight of all paths completely unre- 
stricted. Here S = \ Jq Jq 4>( t i)G(ti — T 2 )4>(T2)dTidT2 
is the "action" with G(t\ — T2) being the inverse matrix 
of the Gaussian correlator f(r\ — t 2 ). Since, Pq(T) is 
expected to decay as exp(— 9T) for large T, we get, 



lim i log P (T). 

1 — *oo 1 



(12) 



If we now interpret the time T as inverse temperature (3, 
then 9 = Ei — Eq where E\ and Eq are respectively the 
ground states of two "quantum" problems, one with a 
"hard" wall at the origin and the other without the wall. 

For concreteness, first consider the Markov process, 
f(T) = exp(— A|T|). In this case, it is easy to see that S 
is the action of a harmonic oscillator with frequency A. 
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The ground state energy, Eq = A/2 for an unrestricted 
oscillator with frequency A. Whereas, for an oscillator 
with a "hard" wall at the origin, it is well known that 
Ei = 3A/2. This then reproduces the Markovian result, 
6 = Ei — Eq = A. For processes close to Markov process, 
such that f{T) = exp(— AT) + e/i(T), where e is small, 
it is then straightforward to carry out a perturbation ex- 
pansion around the harmonic oscillator action in orders 
of e Q. The exponent 6, to order e, can be expressed as, 

9\ f°° 

e = X (l-e— A(T)[l-exp(-2AT)]- 3 / 2 dT). (13) 
n Jo 

At this point, we go back momentarily to the zero tem- 
perature Glauber dynamics of Ising model. Note that the 
spin at a site in the Ising model takes values either 1 or 
— 1 at any given time. Therefore, one really cannot con- 
sider the single site process s(t) as a Gaussian process. 
However one can make a useful approximation in order 
to make contact with the Gaussian processes discussed 
so far. This is achieved by the so called "Gaussian clo- 
sure" approximation, first used by Mazenko in the 
context of phase ordering kinetics. The idea is to write, 
s(t) = sgn[4>{t)] where <p(t) now is assumed to be Gaus- 
sian. This is clearly an approximation. However, for 
phase ordering kinetics with nonconserved order param- 
eter, this approximation has been quite accurate [ p9| . 
Note that, within this approximation, the persistence or 
no flipping probability of the Ising spin s(t) is same as the 
no zero crossing probability of the underlying Gaussian 
process 4>(t). Assuming cj)(t) to be Gaussian process, one 
can compute its two-point non-stationary correlator self- 
consistently. Then, using the same "log-time" transfor- 
mation (with T = log(i)) mentioned earlier, one can eval- 
uate the corresponding stationary correlator f(T). We 
are thus back to the general problem of zero crossing of a 
GSP even for the Ising case, though only approximately. 

In 1 dimension, the correlator f(T) of the un- 
derlying process can be computed exactly, f(T) = 
y/2/(l + exp(2|T|) H and in higher dimensions, it can 
be obtained numerically as the solution of a closed dif- 
ferential equation. By expanding around, T = 0, we find 
that in all dimensions, a = 1 and hence they represent 
non-smooth processes with infinite density of zero cross- 
ings. Hence we can not use IIA or series expansion result 
for 9. Also due to the lack of a small parameter, we can 
not think of this process as "close" to a Markov process 
and hence can not use the perturbation result. However, 
since 9 = Ei — Eq quite generally and since a — 1, we can 
use a variational approximation to estimate Ei and Eq. 
We use as trial Hamiltonian that of a harmonic oscillator 
whose frequency A is our tunable variational parameter 
j| . We just mention the results here, the details can be 
found in For example, in d = 1, we find 9 w 0.35 
as compared to the exact result 9 = 3/8. In e? = 2 and 
3, we find 9 « 0.195 and 0.156. The exponent in 2-d 
has recently been measured experimentally H in a liquid 



crystal system which has an effective Glauber dynamics 
and is in good agreement with our variational prediction. 

So far we have been discussing about the persistence 
of a single spin in the Ising model. This can be imme- 
diately generalized to the persistence of "global" order 
parameter in the Ising model [ljj ]. For example, what is 
the probability that the total magnetization (sum of all 
the spins) does not change sign upto time t in the Ising 
model? It turns out that when quenched to zero temper- 
ature, this probability also decays as a power law ~ <~ 09 
with an exponent 9 g that is different from the single spin 
persistence exponent 9. For example, in 1-d, 9 g = 1/4 
exactly @ as opposed to 6 = 3/8 j|. A natural in- 
terpolation between the local and global persistence can 
be established via introducing the idea of "block" per- 
sistence [Q. The "block" persistence is the probability 
Pi(t) that a block of size I does not flip upto time t. As 
I increases from to oo, the exponent crosses over from 
its "local" value 9 to its "global" value 9 g . 

When quenched to the critical temperature T c of the 
Ising model, the local persistence decays exponentially 
with time due to the flips induced by thermal fluctuations 
but the "global" persistence still decays algebraically, 
~ t~ 8c where the exponent 9 C is a new non-equilibrium 
critical exponent |ll| . It has been computed in mean field 
theory, in the n — > oo limit of the 0(n) model, to first 
order in e = 4 — d expansion p"l[ . Recently this epsilon 
expansion has been carried out to order e 2 |l2| . 

Recently the persistence of a single spin has also been 
generalized to persistence of "patterns" in the zero tem- 
perature dynamics of 1-d Ising or more generally g-state 
Potts model. For example, the survival probability of 
a given "domain" was found to decay algebraically in 
time as ~ t~ B d fl4|| where the q-dependent exponent 
9 d (2) « 0.126 U| for q = 2 (Ising case), different from 
6 = 3/8 and 6 = 1/4. Also the probability that a "do- 
main" wall has not encountered any other domain wall 
upto time t was found to decay as ~ t~ 01 with yet another 
new exponent 6i{q) where 0i(2) = 1/2 and #i(3) ps 0.72 
p5| . Thus it seems that there is a whole hierarchy of non- 
trivial exponents associated with the decay of persistence 
of different patterns in phase ordering systems. 

Another direction of generalization has been to investi- 
gate the "residence time" distribution, whose limiting be- 
haviour determines the persistence exponent [ [jof . Con- 
sider the effective single site stochastic process cj>(t) dis- 
cussed in this paper. Let r(t) denote the fraction of time 
the process (f>(t) is positive (or negative) within time win- 
don [0,i\. The distribution f(r,t) of the random vari- 
able r is the residence time distribution. In the limits 
r — ► or r — * 1 , this distribution is proportional to usual 
persistence. However the full function f(r,t) obviously 
gives more detailed information about the process that 
its limiting behaviours. This quantity has been stud- 
ied extensively for diffusion equation poLEnJ , Ising model 
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p2| , Le'vy processes |3J|, interface models |3(| and gen- 
eralized Gaussian Markov processes |$7| . 

The various persistence probabilities in pure systems 
have recently been generalized to systems with disorder 
[pTsfl . For example, what is the probability that a random 
walker in a random environment (such as in Sinai model) 
does not cross the origin? Analytical predictions for the 
persistence in disordered environment have been made 
recently based on an asymptotically exact renormaliza- 
tion group approach |l8| . 

Another important application of some of these persis- 
tence ideas, experimentally somewhat more relevant per- 
haps, is in the area of interface fluctuations |34|,|35|]. The 
persistence in Gaussian interfaces such as the Edwards- 
Wilkinson model, the problem can again be mapped to 
a general GSP but with a non-Markovian correlator . 
In this case, several upper and lower bounds have been 
obtained analytically BJ] . For nonlinear interfaces of the 
KPZ types, one has to mostly resort to numerical means 
p5| . The study of history dependence via persitence has 
provided some deeper insights in the problems of inter- 
face fluctuations |3^j36|| . 

On the experimental side, the persistence exponent has 
been measured in systems with breath figures soap 
bubbles |]ltj|] and twisted nematic liquid crystal exhibit- 
ing planar Glauber dynamics 0] . It has also been noted 
recently [ j38f that persistence exponent for diffusion equa- 
tion may possibly mbe measured in dense spin-polarized 
noble gases (Helium-3 and Xenon-129) using NMR spec- 
troscopy and imaging [|39| . In these systems the polariza- 
tion acts like a diffusing field. With some modifications 
these systems may possibly also be used to measure the 
persistence of "patterns" discussed in this paper. 

In conclusion, persistence is an interesting and chal- 
lenging problem with many applications in the area of 
nonequilibrium statistical physics. Some aspects of the 
problem has been understood recently as reviewed here. 
But there still exist many questions and emerging new 
directions open to more theoretical and experimental ef- 
forts. 
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gestions and discussions. I thank M. Barma, B. Derrida 
and C. Godreche for useful discussions. I am grateful 
to CNRS, Universite' Paul Sabatier for hospitality where 
the whole series of work began. 
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